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Ahstract 

THic naoer describes a method to efficiently and 
accurately approximatetheetfetf of . de ^ 

aJfproStions ^ jjw. mrthJd. 

Differential Equation Based tu / vjbration 

£i» 

perturbing the height, w^h.^cross-sectiona 1 area.^ hp 

mefhod 3 provedto' be very aerate, and in most cases 
was more accurate than the linear iay'ot 

-Tsea ,o 

approximate bending stresses. 
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It is highly desirable in optimization to be able to 

s ->2^ sS£S«j ,in ™ s 

f h e n S ra p?U“^ many opllmization^ procedures 

full analysis can be eliminated. 1 Taylor series and o 
full ana V s ® , s0 used outside of optimization 

tor^aufclT assessn^nt ^fthe S effect of design changes. 

KtcSes and 

are described m r ®* er ®"^3- 1 ^slgn if leant ly 

efficient method is d ^f*®5 it W of approximations 

STa^ihe onSemlal Equation Based (DEB) 
approximations. 

of PER Method 

in this section of the paper the DEB method is 
developed for approximating the effect of de g 
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res P° n se. First, the method will be 
mrtw 9 or approximation of frequencies then 
mode shapes, and finally for static displacements.’ 

Appr oximation Of Freniier^ 

From reference 12, the equation for the derivative 
lo^ a ^ignSe!' b v r s° n elge " value ’ w#h 


J dv dv dv 


and 

D = -<t> T MQ - 0 

2 dv 


(7) 


( 8 ) 


dv L dv dv J 


(1) 


where K is the stiffness matrix, M is the mass matrix and 

is the eigenvector normalized such that <t» T M<& = l 
Define two scalar quantities a and b as follows 


a = <» T ^«*» and b = <t> T ~<t> 

then Equation (1) may be written as 
do) 2 


( 2 ) 


IIS' rnU, » nnrr^ >» '«"* * <*eWW 1 3 . 

equation (6) is a nonlinear first order vector differential 
3£*n in * with variable coefficients. In ofderto 

nomfriafiH equat, ° n ' Q and D evaluated at the 
nominal design and assumed constant. Equations ( 1 ) 

and (5) may be used for explicit forms of inland or in 
the ^calculations of Q and 0. The solutioMo equalion 


4> = Ae Dv -— Q 


(9) 


dv 


■ = b-ao) 2 


(3) 


Equation (3) may be interpreted as a first order 

f‘ ff . er ® n l ial equation in co 2 with variable coefficients If a 
and b do not vary excessively with v, they may be 

constarft^ Thin e th° mmal de , Si9n and conside red to be 

Wairsaor" so,u,ion ,o equa,ion <»>• 


where A is determined from the nominal condition that 

"5® n . v * v Solving for A and substituting into 
equation (9) g,ves the explicit approximation for <6 as 


«=(®o+iQ]e 0 ( , - v ")-iQ 


( 10 ) 


<B 2 = Ce~ av + — 
a 


(4) 


Equation (10) is a vector equation, however it k 
uncoupled in the sense that each component o i 
varies independently with the design variable v The 

^l 6 " C l betWeen com Ponents is reflected In the 
corresponding components of and Q. 


Solving for C from the condition of co 2 = co 2 when v-v 
and substituting into Equation (4) gives an expiicS 
approximation for w 2 as 


Approximation of Qjspiacgmgiiis 


a 2 =f(o%~— )e~ a ^ v-v ° ) + 

\ aj a 

Approximation of Mode 


(5) 


__ Assuming a linear elastic structure in which the 

ssssa saa r wiih * ma s 


dU 

dv =BU 


Of ,o F :°^H refe K renC6 13, the e P uation for the derivative 
of the mode shape, 4 > with respect to a design variable 


where 


d<t* „ ^ 
— = Q + D<t> 
dv 


B = -K~‘ — 
dv 


(11) 


( 12 ) 


( 6 ) 


where Q satisfies the equation 


Equation (11) is a differential equation that can be 
solved explicitly for an approximation of the 
displacement, U. Similar to solving the differential 
txtuanons tor frequencies and mode sha£s 
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B is evaluated at the nominal design and assumed 
constant for moderate excursions from the nominal 
design. Equation (11) is a vector differential equation 
having the following solution 

U = e Bv A ( 13 ) 


where e B is the matrix exponential defined as 


b B 2 B 3 

e 3l+B+ 3r + ir 


(14) 


Solving for A in equation (13) from the condition that 
U=U 0 at v=v 0 gives 


Also pointed out in reference 14, equations (16) and 
(19) differ only by the factorial in each term. Equation 
(16) is less accurate than the Taylor series equation (19) 
when they both converge. Equation (16) is based on 
the power series for the exponential and will converge 
for any value of Av, whereas the Taylor series may not 
always converge. Since the terms in both equations 
(16) and (19) are easily evaluated using equation (17), 
the DEB method would not be the preferred 
approximation unless the Taylor series did not 


U = e B < v - v °)u 0 < 15 ) 

Using equation (14) in equation (15) gives the 
displacement approximation 

u=u„*B(v-v 0 )u (16) 

Rather than compute the matrix B and its powers and 
then multiply by U 0 , the B n U 0 terms were obtained 
iteratively as follows 


converge. 

DEB Method Extended to Multiple Design Variables 

The approximations in equations 5,10, and 15 are 
directly applicable only when a single design variable is 
perturbed. To calculate the effect of simultaneous 
changes in several design variables, the DEB method 

is modified as follows. If v o is the vector of nominal 
design variables and v is the vector of perturbed 

design variables, then v-v o is the direction of the 
move from the nominal to the perturbed design. Any 
design along the move direction may be associated with 
a scalar 0 between 0 and 1 . Then 


KBU 0 

KB 2 U 


dK 

dv 


U„ 


dv 


KB n U 0 = -^B n_I U 0 


The Taylor series expansion of U is 


v(0) = v o + 0(v-v o ) 


( 20 ) 


(17) In particular, v(l) = v and v(0) = v 0 . Now the sensitivity 

equations (1). (6). and (11) may each be written with 9 
as the independent variable. For example, equation (1) 
is rewritten as 



* T 




or 


(21a) 


dU, , li^U, N 2 
u = ll<+ _ (v . Vo)t __ r (v-v 0 )*. 


(18) 


dor . , 2 

_ = b e -a 9 a> 


(2ib) 


As pointed out in reference 14, when K is a linear 


where 


function of v, equation (18) can be written as 
U = U 0 + BU 0 (v - v 0 ) + B 2 U 0 (v- v 0 ) 2 +- 


(19) 


*TdM . 

a A = <t> — <& 
6 d0 


and b fl =«I> T — <fr 
u d0 


( 22 ) 


and 
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(23) 


v(e,y) = V 0 + e(v- v„) = v(y o )+0[v(y)-v(y # )] (27) 


dM " dM dvj 
d0 ~ i,i dvj d0 
dK = " dKdvi 
d0 iTidVj d0 




where n is the number of design variables. Now the 
solution to equation (21 b) is 


a, 2 (0) = 


* 2 be ' ' 

(On *■ 

a e ) 


} -a 9 (8-8 0 ) + bg_ 
a 8 


(24) 


The method for nonlinear design variables then 
becomes very similar to the DEB method for multiple 
linear variables (equation 26) . For this case 


dM 

d0 

dK 

d0 


" dM dvj 
iTidVj d9 
"dKd^ 

i-idvj d0 



(28) 


Using equation (20) co 2 (l) = co 2 (v) and by definition 

0 o = O and 0=1 at the perturbed design. Therefore, 
the final form of equation (24) is 

(o 2 (v) = f(o 2 -^V a9+ — (25) 

l a 9 J a 9 

Equation (25) corresponds to one step |v — v 0 | . For a 
step size of sjv - v 0 | 

(0 2 ( v ) = ft0o e -aeS + -^2- (26) 

v a e ; a e 


Then a 9 and b 9 are calculated from equation (22). 

Unlike the method for linear variables, a 9 and b 9 are 
evaluated at each step along the move direction due to 
the nonlinearity of the perturbed variable and in 

equation (26) s is always equal to 1 . Evaluating a 9 and 

b 9 involves a recalculation ofvj(y) for each perturbed 
step of y and the multiplication of a matrix by a scalar. 
The derivatives of the mass and stiffness matrices do 
not have to be recalculated. They are evaluated on c e at 
the nominal design and used throughout the range of 
the move direction. Again, a similar extension has been 
carried out for and U. 


Extension of the method to multiple design variables 
has also been carried out for <6 and U. 

DEB Method Extended To Nonlinear Design Variables 

We will define a nonlinear design variable to be one 
that contributes nonlinearly to either the mass or 
stiffness matrix and a linear variable is one that 
contributes linearly to both the mass and stiffness 
matrices. An early attempt at applying the DEB method 
in a straight forward manner to nonlinear design 
variables proved to be highly inaccurate. In retrospect, 
the reason for the inaccuracy was that the coefficients 
which were assumed constant (e.g. a and b) were in fact 
highly nonlinear functions of the design variable. The 
solution to this difficulty is as follows. Define y to be a 
vector of nonlinear design variables and v to be a 
vector of linear variables. It is assumed that the linear 
variables may be written as functions of the nonlinear 

design variable such that v 0 = v(y 0 )and v = v(y). For 
example, the bending inertia I and the cross-sectional 
area A contribute linearly to M and K and they are both 
functions of the nonlinear variable H. Then equation 
(20) may be rewritten with v as a function of y and 0 


Results 

The approximations for frequencies, mode shapes, 
and displacements were tested on a 193 inch 
cantilever beam modeled with ten equal elements (see 
fig.l). The beam has a box cross-section and a tip 
mass. Properties of the beam are listed in table 1 . 




Fig.l Test problem. 
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Table 1 Material properties and cross-sectional 
dimensions of beam model 


Element No. 

1 

2-10 

E (psi) 

4.90E6 

5.85E6 

I P(lb/in 3 ) 

0.07 

0.07 

H (in.) 

5.00 

5.00 

B (in.) 

3.75 

3.75 

t (in.) 

0.80 

0.80 

d (in.) 

0.10 

0.10 


The test cases involved perturbing the height H, width 
B, cross-sectional area A, tip mass M, and bending 
inertia I of the beam. All calculations were generated in 
a general purpose finite-element program (EAL) 15 . 
Results were compared with exact solutions and the 
linear Taylor series approximation. 

Frequency A pproximation Results 

The DEB method shows excellent capability for 
approximating frequencies of perturbed designs. For 
example, figure 2 shows a graph of the fundamental 
eigenvalue versus the height of the beam, H for values 
perturbed from the nominal value of 5.0 inches. For as 
much as a 50 percent increase in H, the new 
approximation is within 2 percent of the exact solution 
compared to 12 percent for the linear Taylor series 
approximation. It is also evident from the figure that 
beyond a 45 percent decrease in H, the linear Taylor 
series method gives negative values for the 
eigenvalues. This is an example of the DEB method 
applied to a nonlinear variable. Figure 3 shows a 


Linear Taylor series 



Fig. 2 Comparison of DEB and Taylor series 
approximations of frequency for perturbation of beam 
depth, H ( B-3.75 in., M=»0. Ibm. ). 

comparison between the approximations for the 
frequency due to a change in the magnitude of the tip 
mass. Again the new approximation is very accurate. 


Only a four percent deviation from the exact solution is 
seen for a 100 percent increase in the tip mass 
compared to a ten percent deviation for the linear Taylor 
series approximation. Figure 4 illustrates the application 



Fig. 3 Comparison of DEB and Taylor series 
approximations of frequency for perturbation of tip mass 
( B=3.75 in., H=5.0 in. ). 

of the DEB method to simultaneous changes in three 
design variables: tip mass, bending inertia, and cross- 
sectional area of the beam. Each step, s represents a 
ten percent perturbation in the design variables from 
their nominal values. In comparing the DEB 
approximation with the exact solution there was a 36 
percent error for 50 percent increase in the design 
variables. The linear Taylor series approximation was in 
error by nearly 100 percent. Figure 5 shows the DEB 
and linear Taylor series approximations of the second 
bending frequency compared to the exact solution for 
perturbed values of the beam depth, H. Again, the new 
method is very accurate. For as much as a fifty percent 
increase in H, the DEB approximation is within 2 percent 
of the exact solution compared to 12 percent for the 
linear Taylor series. This result is very similar to the result 
for the first bending mode with H as the design variable. 
Similar trends occurred in perturbation studies of the 
other variables for approximating frequencies of the first 
and second bending modes. 
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Number of steps, s 

Fig. 4 Comparison of DEB and Taylor series 
approximations of frequency for simultaneous 
perturbation of tip mass M, bending inertia I, and cross- 
sectional area A ( Nominal values: M-5.0 Ibm., 1-28 863 
in 4 , A-7.0 in 2 ). 


Figure 7 shows the second bending mode shape for 
the nominal design and both approximations and the 
exact shape for an 80 percent increase in the tip mass. 
Although the linear Taylor series approximation is very 
accurate, the DEB approximation is slightly closer to the 
exact curve in both figures. However, the mode shapes 
do not seem to be sensitive to design changes for this 
test problem as indicated in figures 6 and 7. 



Normalized distance along beam 



Fig. 5 Comparison of DEB and Taylor series 
approximations of frequency for perturbation of beam 
depth, H (B-3.75 in., M-0. Ibm. ). 

Mode Shape Approximation Results 

Figure 6 shows the first bending mode shape for 
the nominal design, and the linear Taylor series 
approximation, the DEB approximation, and the exact 
solution for a 50 percent increase in the beam width, B. 


Fig. 6 Comparison of DEB and Taylor series 
approximations of first bending mode shape for 50 
percent increase in beam width, B ( Nominal values: 
H-5.0 in., B-3.75 in. ). 



Normalized distance along beam 


Linear Taylor Series 
DEB method 
Exact solution 
Nominal 


Fig. 7 Comparison of DEB and Taylor series 
approximations of second bending mode shape for 80 
percent increase of tip mass, M ( Nominal values: 
B-3.75 in., H-5.0 in., M-5.0 Ibm. ). 
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In some cases the DEB method did not approximate the 
mode shapes as accurately as the linear Taylor series 
approach. Usually, these were cases in which the 
change in the mode shape due to perturbations in the 
design were small. Other test cases may provide a 
more challenging problem for testing these 
approximation methods for mode shapes and will be a 
topic for future investigation. 

Static Displacement Ap proximation Results 

The cantilever beam of figure 1 was also used to 
test the method for approximating static displacements. 
The design variables were the dimensions H and B. A 
cubically distributed loading was used to represent the 
applied force and the load values at each grid point of 
the model are given in table 2. 

Table 2 Values of applied force at the grid points of 
the beam model 


Grid point 

Force (Ibf) 

Grid point 

Force (Ibf) 

1 

0.000 

7 

5.400 

2 

0.025 

8 

8.575 

3 

0.200 

9 

12.800 

4 

0.675 

10 

18.225 

5 

1.600 

1 1 

25.000 

6 

3.125 




Figure 8 shows a typical convergence plot for the DEB 
method equation (16) for the maximum displacement. 
This figure suggests that the number of terms needed 
for a converged value is 3 or 4. For this same number of 
terms and the same amount of computation, the Taylor 
series approximation (equation 19) is known to be more 
accurate and would be chosen over the DEB method 
unless it did not converge for certain values of Av. 
Results for the linear Taylor series will be included 
because it is a frequently used method of 
approximation. 

Figures 9, 10, and 11 show results for perturbation 
of the nonlinear design variable, H. In figure 9, the 
maximum displacement for values of H from 3.5 inches 
to 8.5 inches was plotted for the DEB, linear Taylor 
series approximations, converged Taylor series and the 
exact solutions. As H is decreased, the bending 
stiffness of the beam is decreased causing the beam 
displacement to become very large. At this end of the 
range it is difficult for the approximation methods to 
follow the curvature of the exact results. At the other 
end, however, the DEB method appears to follow the 
trend of the exact curve. The linear Taylor series 
approximation gives negative values for the 
displacement beyond 40 percent perturbation of H from 
the nominal value. 



Fig. 8 Convergence plot of DEB method 



Fig. 9 Comparison of DEB and Taylor series 
approximations of maximum displacement versus 
perturbation of beam depth, H. 
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Figure 10 is a graph of the displacement of the beam for 
a 30 percent decrease in the nonlinear variable H. The 
figure compares the DEB and linear Taylor series 
approximations to the exact solution. Both 
approximations are very good with the DEB method 
slightly more accurate. 



Normalized distance along beam 


Fig. 10 Comparison of DEB and Taylor series 
approximations of displacement for 30 percent 
decrease in beam depth, H. 


A similar comparison is shown in figure 11. For a 50 
percent increase in H from the nominal value, the DEB 
method had an error of 47 percent for the maximum 
displacement. However, the linear Taylor series 
approximation gives a displacement in the opposite 
direction. This is in accordance with the result we saw in 
figure 9. 



Normalized distance along beam 


Fig. 1 1 Comparison of DEB and Taylor series 
approximations of displacement for 50 percent increase 
in beam depth, H. 

Figures 12, 13, and 14 show results for linear Taylor 
series and DEB displacement approximations, 
converged Taylor series and exact solutions for the 
linear variable B. In figure 12, the maximum 
displacement for values of B from 1.875 inches to 
5.625 inches was plotted for each approximation 
method and compared to the exact solution. Again, it 
appears to be more difficult to approximate the 
displacement for small beam dimensions, but the DEB 
method follows the exact curve, particularly for the 
larger values of B. Figure 13 shows a plot of the 
approximate displacements and the exact solution for a 
50 percent decrease in B. Both approximations are 
good but the DEB curve is closer to the exact. Figure 
14 shows the result for a 50 percent increase in the 
design variable. It shows that the DEB approximation is 
within 8 percent of the exact solution compared to the 
linear Taylor series error of 21 percent. 
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Displacement, U(in.) Maximum displacement, (in.) 



Fig.12 Comparison of DEB and Taylor series 
approximations of maximum displacement versus 
perturbation of beam width, B. 



Fig. 13 Comparison of DEB and Taylor series 
approximations of displacement for 50 percent 
decrease in beam width, B. 


Linear Taylor series 



Normalized distance along beam 


Fig. 14 Comparison of DEB and Taylor series 
approximations of displacement for 50 percent increase 
in beam width, B. 

Additional Approximations 

An additional application of the DEB approximation is 
in estimating system response quantities that depend 
on approximated quantities. Inserting the approximate 
displacement field into the stress-displacement 
equations to obtain approximate stress field is an 
example of this. Figure 15 shows a graph of maximum 
normal bending stress for perturbed values of B. The 
stresses were obtained using the approximate values of 
U from the DEB method and the linear Taylor series 
approach and are compared to the converged Taylor 
series and exact solutions. The approximate stress field 
shown in figure 15 is very accurate. This was a 
particularly good result in view of the fact the the 
stresses are in effect proportional to the second 
derivatives of the approximate displacements. 
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« \ Linear Taylor series 
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tT 21. - x \ DEB method 
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0 2.5 3.5 4.5 5.5 6.5 

Beam width, B (in.) 



0 1 2 3 4 5 

Number of steps, s 


Fig. 15 Comparison of DEB and Taylor series 
approximations of maximum stress versus perturbation 
of beam width, B. 



Recall that the methods developed in this paper were 
based on the assumption that certain variables were 
assumed constant, for example a and b in equations (5) 
and (26). In order to assess the validity of these 

assumptions, the actual values of the coefficients a 9 

and b e were plotted for values of the tip mass, bending 
inertia, and cross-sectional area perturbed 

simultaneously. Figure 1 6 shows a graph of a e and 
versus the number of 1 0 percent step increases in the 

design variables. The lower curve, shows that ^2. is 

a e 

essentially constant. The upper curve representing the 
coefficient a 9 shows a deviation of up to 19 percent 
from the nomial value when the design variables are 
changed by up to 50 percent. Assuming a e to be 
constant was not detrimental to the approximation 
results. As shown in figure 4, the DEB method was able 
to obtain very good approximations for the full range of 
the pertubed design variables. 


Fig. 1 6 Actual values of coefficients assumed constant 
in DEB frequency approximation for simultaneous 
perturbation of tip mass, bending inertia, and cross- 
sectional area 


Figure 17 shows a graph of the coefficients that appear 
in the mode shape approximation of equation (10). 
The plot of D versus perturbed values of tip mass show 
plus and minus 20 percent deviations in D based on up 
to 100 percent perturbation in the tip mass. The 

second plot shows the vector Q for the nominal tip 

D 

mass and for plus and minus 100 percent perturbations 
from the nominal mass. These plots show that the 
vector varied considerably, however, as shown in figure 
6, this did not seem to have a large effect on the 
approximate mode shapes obtained with the DEB 
method. 
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Fig. 17 Actual values of coefficients assumed constant 
in DEB mode shape approximation for mass 
perturbation 


Concluding Remarks 

This paper described a method that is based on 
interpreting sensitivity equations of structural response 
as differential equations that may be solved for closed 
form approximations. The method was developed for 
approximating vibration frequencies, mode shapes and 
static displacements. The test cases used to evaluate 
the approximations involved perturbations of the 
height, width, cross-sectional area, tip mass, and 
bending inertia of a cantilever box beam. The results 
were compared to the linear Taylor series 
approximations and exact solutions. The results 
showed that the new method is very accurate and in 
most cases more accurate than the linear Taylor series 
approximations. The method is applicable to both linear 
and nonlinear design variables and to simultaneous 
perturbation of several design variables. Also, the 
approximated response may be used to approximate 
behavior derived from the response, for example, 
approximating stresses from an approximate 
displacement field. To date, the method has been 
developed and demonstrated for frequencies, mode 
shapes, static displacements, and stresses derived 
from the displacement. In principle, the method is 
applicable to approximating any quantity for which an 
analytical sensitivity formula is available. 
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